Alignment-Dependent Ionization of Molecular Hydrogen in Intense Laser Fields 
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The alignment dependence of the ionization behavior of H2 exposed to intense ultrashort laser 
pulses is investigated on the basis of solutions of the full time-dependent Schrodinger equation within 
the fixed-nuclei and dipole approximation. The total ionization yields as well as the energy-resolved 
electron spectra have been calculated for a parallel and a perpendicular orientation of the molecular 
axis with respect to the polarization axis of linear polarized laser pulses. For most, but not all 
considered laser peak intensities the parallel aligned molecules are easier to ionize. Furthermore, it 
is shown that the velocity formulation of the strong-field approximation predicts a simple interference 
pattern for the ratio of the energy-resolved electron spectra obtained for the two orientations, but 
this is not confirmed by the full ab initio results. 

PACS numbers: 32.80.Rm, 33.80.Rv 



Time-resolved imaging of the dynamics of nuclei and 
electrons on a femtosecond or even sub-femtosecond time 
scale is a prerequisite for the real-time investigation of the 
formation and breaking of chemical bonds. Ultrashort 
laser pulses have recently been demonstrated to pave a 
possible path to the experimental realization of this long- 
standing dream. For example, ways have been proposed 
and experimentally verified that allow monitoring nuclear 
motion with sub-femtosecond and sub- Angstrom resolu- 
tion in real time [H, [3i It was also experimentally 
demonstrated that the high-harmonic radiation or the 
electrons emitted in an intense laser pulse may in prin- 
ciple reveal information on the electronic structure 0, 0] 
and thus have the potential for time-resolved imaging of 
changes of the electronic structure in, e.g., a chemical 
reaction. To reach this goal it is, however, important 
to understand the relation between electronic structure 
and the strong-field response of molecules. This includes 
the fundamental question whether the rather clear corre- 
spondence between the symmetry of the highest-occupied 
molecular orbital (HOMO) and the strong-field signal as 
indicated for N2 and O2 in 0, Q is really a universal 
phenomenon. 

Already some time ago it was found that within the 
molecular strong-field approximation (MO-SFA) — for- 
mulated in velocity gauge (VG) and within the frame- 
work of a linear combination of atomic orbitals (LCAO) 
— the molecular response to intense laser fields should 
reflect the symmetry of the highest occupied molecular 
orbital (HOMO) In the case of diatomic molecules, 

a simple interference picture arises in the MO-SFA- VG 
that seems to plausibly explain the occurrence or absence 
of suppressed ionization [G]. The term suppressed ioniza- 
tion describes the effect that a molecule with the same 
ionization potential as the one of some so-called compan- 
ion atom is harder to ionize in an intense laser pulse. 
For example, molecular nitrogen shows a similar ioniza- 
tion behavior as atomic Ar, while the ionization yield of 



oxygen is much smaller than the one of Xe atoms; de- 
spite the almost identical ionization potentials of either 
N2 and Ar or O2 and Xe. Although suppressed ioniza- 
tion is also predicted by tunneling models like molecular 
Ammosov-Delone-Krainov (MO-ADK) Q or the length- 
gauge formulation of the MO-SFA these theories do 
not provide a simple interference picture for the effect. 
On the other hand, the energy-resolved electron spectra 
measured in seemed to further support the concept 
of of symmetry-induced quantum-interference effects as 
predicted by MO-SFA- VG. 

One key step towards time-resolved imaging is the 
measurement of the molecular strong-fleld response 
within a molecule-fixed coordinate system. This is also 
the basic difference between the electron spectra mea- 
sured in [5| in comparison to the ones in m. While 
the former are obtained as a function of the alignment 
between the molecular axis and the laser- field axis, the 
latter are averaged over all orientations. This explain the 
increased recent interest even in alignment-dependent to- 
tal ionization yields of molecules in intense laser pulses 
l3, In fact, the results in 11 1 seem to indicate 
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that, at least for N2 and O2, structural information (on 
the HOMO) may be obtained already from such integral, 
but angular-resolved ionyields. Besides the problem that 
the results obtained in |ll| for CO2 seem to be difficult to 
interpret as a simple mapping of its HOMO, already rel- 
atively simple diatomic molecules like O2 and N2 possess 
a HOMO that is not necessarily well described as a linear 
combination of two atomic orbitals, as is easily seen from 
the recent debate about the prediction of MO-SFA- VG 
for the parallel to perpendicular strong-field ionization 
yield of N2 (see ^] and references therein) . Therefore, 
molecular hydrogen with a comparatively simple orbital 
structure and the lack of core orbitals that may disturb 
the strong-field response appears to be a perfect candi- 
date to investigate whether the interference effects pre- 
dicted by MO-SFA- VG occur. Furthermore, H2 has at 
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FIG. 1: (Color online) Ionization yields of H2 (internuclear 
separation R = 1.4 no) for a 10-cycle cos'^-shaped 800 nm 
laser pulse and a parallel (blue circles) or perpendicular (red 
squares) alignment with respect to the laser polarization vec- 
tor. 



least for fixed nuclei now become accessible to in princi- 
ple exact theoretical calculations, even for a non-parallel 
orientation of the molecular axis with respect to the field 
that requires a full 6-dimensional treatment ^ISj . Its ex- 
tension to laser fields with a wavelength of 800 nm and 
the extraction of energy-resolved electron spectra is re- 
ported in this work. 

In fact, the purpose of this Letter is threefold. First, 
the alignment dependence of the total ionization yield 
of H2 in intense laser pulses with a wavelength of 
about 800 nm (Ti:sapphire) is investigated by means 
of a full solution of the time-dependent Schrodinger 
equation (TDSE) describing both correlated electrons in 
full dimensionality within the fixed-nuclei and the non- 
relativistic dipole approximations. It is shown that the 
ratio of ionization yields for parallel and perpendicular 
alignment is a non-monotonic function as a function of 
laser peak intensity, even if focal- volume averaging is per- 
formed, and is in reasonable agreement to a recent exper- 
iment 



12l |. Second, it is demonstrated theoretically that 



MO-SFA-VG predicts a clear interference pattern in the 
ratio of energy-resolved electron spectra for a parallel and 
a perpendicular alignment of a molecule like II2 . Finally, 
the corresponding energy-resolved electron spectra are 
extracted from the TDSE calculation and it is demon- 
strated that the simple interference pattern predicted by 
MO-SFA-VG is not confirmed. 

Figure [1] shows the ionization yields for a parallel and 
a perpendicular orientation of H2 exposed to a 800 nm 
laser pulse as a function of laser peak intensity (for a 
fixed internuclear separation R — 1.4 oq). As in all cal- 
culations shown in this work 10-cycle cos^-shaped laser 
pulses (~ 10 fs FWHM) are used. The numerical ap- 
proach for solving the TDSE of H2 has been described 
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FIG. 2: (Color online) Ratio of parallel to perpendicular ion- 
ization yields of H2 (7i = 1.4ao) as a function of laser peak 
intensity for different wavelengths and 10-cycle cos^-shaped 
laser pulses: without (left panel) and with (right panel) focal- 
volume averaging. For a better comparison to the experimen- 
tal data from 12] the 800 nm result is also shown scaled down 
by a factor 1.18. Furthermore, the (simplified) MO-SFA-LG 
data in [l2| are depicted. 



recently [l5|. It is based on an expansion of the time- 
dependent wavefunction in terms of box-discrctized field- 
free eigenstates obtained from a configuration-interaction 
calculation based on eigenstates calculated in a B- 
spline basis. According to Fig. [T] the ionization yields 
show evident structures from resonantly-enhanced mul- 
tiphoton ionization (REMPI) and channel closings, de- 
spite the rather high intensity and long wavelength. Since 
these structures depend on the orientation (the selection 
rules lead to different REMPI intermediate states and 
the different polarizabilities to different shifts of the mul- 
tiphoton thresholds), the results for parallel and perpen- 
dicular orientation vary differently as a function of inten- 
sity. This is especially evident, if the ratio between par- 
allel and perpendicular alignment is considered (Fig. [2]). 
This ratio shows a highly oscillating structure. Notewor- 
thy, in a certain intensity range the perpendicular align- 
ment leads to a larger ionization yield that the perpen- 
dicular one. A small variation of the wavelength within 
780 and 810 nm confirms that even the fine details of 
the structures are reproducible and not some numeri- 
cal artifact. Although focal- volume averaging damps the 
sharp structures, the ratio remains non-monotonous and 
in some peak intensity regime the perpendicular orienta- 
tion is easier to ionize than the parallel one. 

Very recently, the ratio of the ionization yields for 
parallel and perpendicular alignment was measured for 
800 nm laser pulses in [l3| . The results are also shown in 
Fig. [21 Note, that the intensities have been adopted as 
in T2i , although the experiment was performed with cir- 
cul ar p olarized light. This is in agreement to the analysis 
in [l2| that compared to a model for linear polarization 
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using unsealed intensities. Supported also by the present 
finding, the authors of jU] argued that due to experimen- 
tal reasons the found anisotropy may be underestimated. 
Dividing the TDSE results by a factor 1.18 gives in fact 
a very satisfactory agreement between theory and exper- 
iment. Although not really statistically relevant in view 
of the error bars, one may note that also the experimental 
data indicate a small local maximum whose exact posi- 
tion depends according to the present TDSE results on 
the exact wavelength. In [13] a simphfied MO-SFA-LG 
model is also proposed that considers basically the differ- 
ent polarizabilities of H2 for parallel and perpendicular 
orientation. The model is in the intensity regime shown 
in [1^ in good overall quantitative agreement with the 
present TDSE result, but does not show the small min- 
imum. More importantly, it appears unlikely that the 
model could explain the sharp decrease and even inver- 
sion of the ratio for slightly smaller intensities. 

In the case of a homonuclcar diatomic molecule with 
internuclear separation R a bonding HOMO $ built from 
s-type atomic orbitals </) is given within the LCAO as 



$(r, R) - a{0(r, -R/2) + ^(r, R/2)} 



(1) 



where r is the electronic coordinate and a the normaliza- 
tion constant. According to MO-SFA-VG [1] (also called 
first-order IMST) this leads to the A^-photon ionization 
rates (integrated over the directions of the emitted elec- 
tron) 



= N,, I dk 



(2) 



N 



in a linearly polarized laser field, if the HOMO is occupied 
by Ne electrons. The in Eq. ([2]) occurring differential 
ionization rates are given by 



N) 



dk 



= 2 



N 



X 2a(kAr|(/)) C0S^(kAr • R/2) 



(3) 



Here, — Nuj — {Up + Eion) is the kinetic energy of an 

electron after absorption of N photons, Up = F^/(4w^) is 
the ponderomotive energy of an electron in a laser field of 
frequency cu and peak field strength F, Eion — /2\s the 
ionization energy of the molecule, and = (k^/F)^/*^ is 
a Coulomb correction factor. Finally, (k|(/)) is the Fourier 
transform of the atomic orbital 0(r) and J„(a, 6) is a 
generalized Bessel function of two arguments as defined 
in [1^ . The polarization axis e enters the ionization rate 
only through one of the arguments of the Bessel function, 
ao = (F/a;2)e. 

The sum of the Af-photon ionization rates for all en- 
ergetically allowed values of A^ yield the energy-resolved 
electron spectra also known as above-threshold ionization 
(ATI) spectra. Since the Fourier transform of spherically 



symmetric (s-type) orbitals is also spherically symmet- 
ric, it depends only on k. Therefore, the ratio of the 
A'^-photon ionization rates for parallel and perpendicular 
orientations of the field polarization vector with respect 
to the internuclear axis can be written as 



Xjsi — 



pl^ / dk ^.gAre • k, 6^ cos^ ^dAfk ■ 



R 



^ N 



J dk.J'^(^gNe ■ k, bj cos^^dArk ■ 



R 



(4) 



where — aokj^, b — Up/{2uj), d^ — Rkj^/2, and all 
factors depending only on the absolute value of k^r were 
taken out of the integral and cancel each other when the 
ratio is considered. 

Fixing the coordinate system in such a way that 
its z axis agrees with the polarization vector e, one 

has r" = (0,0,1), R^ = (1,0,0), and k = 
(sin 6* cos 0, sin sin (/), cos 0). Then the ratio ^ can be 
rewritten as 



X 



N 



dO sin Jff(^gpj cos 9, cos^ {dN cos 9) 

gM cos 6', 6) [1 + Jo(2dAr sin6')]/2 

(5) 



where the identity 



/ d0cos2(5 COS0) = 7r[l -I- M25)] 
Jo 



(6) 



for zero-order Bessel function of the first kind was used. 

For very strong fields the function jfj{gN cos9,b) 
peaks usually sharply around cos 6* = ±1 (i.e. in the 
case of ionization from a spherically symmetric atom 
the vast majority of electrons are ejected parallel to the 
laser polarization axis). On the other hand, the func- 
tions cos^(dAr cos 9) and Jo{2dN sin^) vary rather slowly 
in these regions. They may thus be approximated by 
their values at cos^ ±1: 



COS^(dAr cos6') « cos^(d 



N 



Jo(2dNsm9) w 1 . (7) 



Substitution of Eq. ([7]) into Eq. ([5]) results finally in a 
very simple expression for Ajv, 



Xn ~ C0S^(d7v) = C0S^(i?fc7v/2) 



(8) 



Note, that for very small d^r (which means small R or 
small k]^) both functions in Eq. ([7]) are equal to 1, so the 
ratio Xn is also 1. 

According to Eq. ^ MO-SFA-VG predicts a pro- 
nounced minimum to occur at the electron energy E « 
7r^/(2i?^), if the ratio between the the energy-resolved 
electron spectra obtained for a parallel and a perpendic- 
ular orientation is considered. Its origin is a pure inter- 
ference phenomenon that is due to the destructive inter- 
ference caused by the phase shift between two electronic 
wavepackets emerging from the two nuclei and moving 
along the polarization axis. 
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FIG. 3: (Color online) Energy- resolved electron spectra for 
parallel and perpendicular orientations of H2 {R = 3.0 ao) 
for a 10-cycle cos^-shaped 800 nm laser pulse with a peak 
intensity of 2 ■ lO" W/cm^. 
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FIG. 4: (Color online) Ratio of parallel to perpendicular 
energy-resolved electron spectra of H2 {R — 3.0 ao) for 10- 
cycle cos^ -shaped 800 nm laser pulse and three different peak 
intensities. The ratios predicted by MO-SFA-VG in Eq.[5]are 
also shown using circles. 



In order to investigate the occurrence of this interfer- 
ence phenomenon predicted by MO-SFA-VG, the energy 
resolved electron spectra of H2 were also extracted from 
the TDSE calculations. In order to amplify the effect, 
the larger internuclear separation R = 3.0 ao is consid- 
ered in Fig. [3l This should be an ideal case, since at 
these distances the HOMO of H2 is extremely well de- 
scribed by two hydrogenic Is orbitals. There are, in 
fact two more reasons for choosing a larger R value. If 
the predicted minimum lies at too high energies, rescat- 
tering that is not incorporated in the first-order SFA 
theory could dominate and cover the interference phe- 
nomenon. Furthermore, the energy-resolved spectra ob- 
tained from the TDSE calculation converge more easily 
for lower energies, due to their fast exponential decay. 



Despite all these attempts to find perfect conditions for 
finding the interference minimum predicted by MO-SFA- 
VG, the ratio of the electron spectra for the two orien- 
tations does not show the behavior predicted by MO- 
SFA-VG in Eq. ^ (that agrees well to the approximate 
expression ([H])) as can be seen from Fig. U) It is inter- 
esting that the TDSE results indicate in fact a rather 
universal overall behavior of the ratio independent of the 
laser intensity. However, the ratio first increases with in- 
creasing electron energy and remains on average above 
unity. As is clear from Fig. [31 the total ionization yield 
stems dominantly from the very low-energy electrons. In 
this regime the ratio between the electron spectra for 
parallel and perpendicular is, however, rather intensity 
dependent. This explains the already discussed inten- 
sity dependence of the total yield found already at the 
equilibrium distance. 
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